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ABSTRACT 

We report on new imaging observations of the Lyman alpha emission line (Lya) , performed with the 
Hubble Space Telescope, that comprise the backbone of the Lyman alpha Reference Sample (LARS). 
We present images of 14 starburst galaxies at redshifts 0.028 < z < 0.18 in continuum-subtracted 
Lya, Ha, and the far ultraviolet continuum. We show that Lya is emitted on scales that systemati- 
cally exceed those of the massive stellar population and recombination nebulae: as measured by the 
Petrosian 20 percent radius, i?p2o, Lya radii are larger than those of Ha by factors ranging from 1 to 
3.6, with an average of 2.4. The average ratio of Lya-to-FUV radii is 2.9. This suggests that much 
of the Lya light is pushed to large radii by resonance scattering. Defining the Relative Petrosian 
Extension of Lya compared to Ha, ^Lya^ ^P20 /^P20' we nn( i ?Lya to be uncorrected with total Lya 
luminosity However ^Lya is strongly correlated with quantities that scale with dust content, in the 
sense that a low dust abundance is a necessary requirement (although not the only one) in order to 
spread Lya photons throughout the interstellar medium and drive a large extended Lya halo. 
Subject headings: Physical data and processes: Radiative transfer — Galaxies: evolution — Galaxies: 
formation — Galaxies: starburst — Cosmology: observations 



1. INTRODUCTION 

The Lyman alpha emission line (Lya), emitted by the 
spontaneous de-excitation over the n = 2 — > 1 elec- 
tronic transition in neutral hydrogen (Hi), is now an 
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established observatio nal probe of evolv i ng galaxies in 
the h igh-z Universe (jCowie JcHxi U99l IRhoads et all 
l2000f ) . Exploita tion of Lya has r e sulted in sign i ficant 
galaxy surveys (lOuchi et all 120081: iNilsson et all 120091: 



IGuaita et al. l l2010HAdams et alJl2011ll the next gener 
ations of which will recover vast numbers of galaxies. 
However the Hi abundance in most galaxies, combined 
with the large Lya absorption cross section of ground- 
state hydrogen, suggests that most Lya will be absorbed 
and re-scattered by the same transition that created it. 
Thus most Lya photons are thought to be subject to 
multiple scattering events as they encounter neu tral gas, 
result ing in a complica ted r adiative transport (iNeufeldl 
fl990t IVerhamme et al.ll20M iLaursen et al.ll2009f ). 

Because Hi is often found at distances that ex- 
ceed the size of stellar disks and star-forming re- 
gions (lYun et allll994tlMeurer et al|[l99ll Cannon et all 
l2004h . characteristic Lya scale lengths may be ex- 
pected to be substantially larger than those of, for ex- 
ample, the FUV continuum or Ha. Indeed th i s has 
been well observed at high z (e.g. Fvnbo et all 120011 : 



Rauch et all 120081; ISteidel et all 12011 



although see also 



Feldmeier et aUl2013D and low ( Ma.s-H csseet all 120031: 



Ostlin et all 120091) . and studied . 
ulation (lLaursen et all 120091: iBarnes fc Haehnel 



extensively by sim- 



iZheng et alll201 hi IVerhamme et alJl2012j ). 

In this Letter we present images from the Lyman alpha 
Reference Sample (LARS). The LARS program (Ostlin 
et al, in prep; Hayes et al, in prep) is targeting 14 UV- 
selected star-forming galaxies in the nearby Universe, all 
of which have been imaged in Lya, Ha, H/3, and five 
UV/optical continuum bands. Many other observations, 
both in hand and ongoing, are providing gas covering 
fractions and kinematics, and measuring the Hi mass and 
extent directly. HST imaging allows us to probe spatial 
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Figure 1. False-color images of the LARS galaxies 01 to 08. Red encodes continuum-subtracted Ha, green the FUV continuum, and blue 
shows continuum-subtracted Lya. Images have been adaptively filtered to show detail. Scales in kpc are given on the side. Intensity scales 
are logarithmic, with intensity cut levels set to show detail. 



scales down to 28 pc in individual galaxies, quantify the 
extent of Lya, and compare it with other wavelengths 
and derived properties. This Letter discusses the exten- 
sion of Lya radiation. In Section [5] we briefly summa- 
rize the data and show the new images. In Section [3] 
we quantify the sizes of the galaxies in Lya, FUV, and 
Ha, and discuss them with reference to high- 2: measure- 
ments in Section [4] In Section [5] we show how a low dust 
content seems to be a necessary prerequisite in order to 
produce this extended emission. We assume a cosmology 
of (H , n u , 0\) = (70 km s" 1 Mpc _x ,0.3, 0.7). 

2. LARS IMAGES 

LARS consists of 14 star-forming galaxies selected by 
FUV luminosity from the GALEX all-sky surveys, and 
imaged with Hubble Space Telescope cameras ACS/SBC, 
ACS/WFC, and WFC3/UVIS. The sample selection, ob- 
servations, and data processing are described in detail in 
Ostlin et al. (in prep). FUV luminosities range between 
log(L FU v/i0) = 9.2 and 10.7, overlapping much of the 
luminosity range of Lyman-brcak Galaxy (LBG) surveys, 
and are listed in Table [1] 

We use the Ly man alpha extraction software (LaXs, 
I Haves et al.ll2009l) to produce continuum-subtracted Lya 
and Ha images, corrected for underlying stellar absorp- 
tion and contamination from [Nil]. In 1 arcsec square 
boxes away from the targets we measure r.m.s. back- 
ground noise of 5.7 x 10~ 19 erg s _1 cm -2 in Lya, 
2.1 x 10~ 21 erg s" 1 cm" 2 A" 1 in the FUV, and 6.8 x 
10~ 19 erg s _1 cm~ 2 in Ha. Total Lya luminosities range 
from (non-detection) and 2 x 10 erg s _1 with a me- 
dian of 8.1 x 10 41 erg s _1 ; roughly seven of the objects 
would be recovered by the deepest Lya surveys (Hayes 
et al. in prep). 

We present our first imaging results in this paper as a 
series of RGB composite images in Figures [T] and [5] In 
green we encode the far UV continuum, which traces the 
unobscured massive stars, and roughly incorporates the 



sites that produce the ionizing photons. In the red we 
show continuum-subtracted Ha, which traces the nebu- 
lae where the aforementioned ionizing photons arc repro- 
cessed into the recombination line spectrum. The con- 
tinuum subtracted Lya observation is encoded in blue. 
The images have been adaptively smoothed using a vari- 
able Gaussian kernel (FILTER/ ADAPTIVE in ES0/MIDAS), 
in order to enhance positive regions of low surface bright- 
ness emission. The intensity scaling of all the images is 
logarithmic, and the levels are set to show the maximum 
of structure and the level at which the faintest features 
fade into the background. 

Immediately it can be seen that Lya morphologies bear 
limited resemblance to those of the FUV and Ha. In 
some cases Lya appears to be almost completely absent: 
LARS 04 and 06 in particular show only small hints of 
Lya emission that contribute negligibly towards filling 
in the global absorption, and the composites arc domi- 
nated by UV and Ha light. Lya is strongly absorbed, 
particularly in the central regions of these objects. Oth- 
ers show copious Lya emission and reveal morphological 
structures that are not seen at other wavelengths. Most 
obviously, LARS 01, 02, 05, 07, 12, and 14 show large- 
scale halos of Lya emission that completely encompass 
the star-forming regions, although the same phenomenon 
is visible to some extent in all the objects, even the ab- 
sorbers. 

We have di scussed this extend e d Lya emission in depth 
in the past (|Haves et al.l 120051 l2007t lAtek et al.l 120081: 
lOstlin et al]|2009D . However now, with an observational 
setup that is more sensitive to faint levels of Lya and a 
larger and UV-selected sample (Ostlin et al, in prep), we 
are able to robustly quantify and contrast these sizes and 
the relative extension of Lya. 

3. APERTURES, SIZES AND GLOBAL QUANTITIES 

In order to quantify the sizes of the galaxies at various 
wavelengths, we adopt the Pctrosian radius (Petrosian 
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Figure 2. Same as Figure [T] except for LARS galaxies 09 to 14. The 

Il976h with index of rj = 0.2: i.e. the radius, R, at which 
the local surface brightness is 20 percent the average sur- 
face brightness inside of R. In Hayes et al. (in prep) we 
will show the Lya extent of some objects to be so large 
that ACS/SBC cannot capture the full flux, and hence 
measurements like 50 percent light radius are not ro- 
bust. Indeed Petrosian radii were developed to be depth- 
independent measures of size. We note from experimen- 
tation, however, that very similar conclusions are reached 
using other definitions. The choice of r\ = 0.2 gives a 
size for every Lya-emitting galaxy in the sample except 
LARS 09, for which even at the full extent of the SBC 
we do not come close to crossing the -q = 0.2 threshold. 
We reach the edge of the detector at 77 ~ 1 (R > 12 kpc) 
and can expect the true extent of Lya to be much larger. 
For the 11 galaxies in which i?p20 ^ s we ^ measured, its 
determination is robust, and would not change were the 
observations deeper or the ficld-of-view larger. i?p2o ; is 
computed for Lya, Ha, and the FUV continuum, and 
listed in Table [TJ Based upon aperture- matched Ha and 
H/3 imaging and standard Case B assumptions, we re- 
cover up to 60 % of the intrinsic Lya flux, although the 
median value is just ~ 3 % (Hayes et al. in prep). 
We compare the light radii graphically in Figure [3l 

The plots show i?p2Q vs - Rp20 ' a comparison that could 
be made at high-z, and Rp2o vs - ^P20> a comparison 
that more directly conveys the difference between the 
observed and intrinsic Lya sizes. Clearly, though, there 
is little difference in the result: Lya radii are, on av- 
erage, substantially larger than corresponding FUV or 
Ha radii. In Table [1] we also report the Relative Pet- 
rosian Extension of Lya compared to Ha, ^Lya, which 

is simply defined as ^P2o/-^P20- galaxies show net 
emission of Lya, where all except for one (LARS 03) 
has ^Lya> 1- The galaxy with the largest extension is 
LARS 14, for which we measure ^Lya=3.6. It is not clear 
whether globally absorbing galaxies LARS 04 and 06 be- 
come emitters on larger scales, but if so their Rp^ must 



black square in LARS 09 masks a UV-bright field star. 

be larger than the radius of the SBC chip, implying that 
^Lya must exceed 5.3 and 13.4, respectively. That would 
make them the most extended objects in the sample. Ex- 
cluding these two galaxies, and also LARS 09 for which 
we can only provide a lower limit, the sample mean (me- 
dian) is computed as 2.43 (2.28). 

4. RELEVANCE FOR HIGH-REDSHIFT STUDIES 

It is important to note that FUV radii imply that all 
the galaxies would be effectively unresolved by ground- 
based observations if they were at z > 2. The largest 
is 8 kpc, which corresponds to the 1 arcsec resolution 
that could be expected from the seeing. However one of 
the objects has a Lya radius of 15.5 kpc: recovering this 
total flux at z ~ 2 would require an aperture of at least 2 
arcsec. Some objects are also highly elongated and were 
they pushed to the high-z Universe, much of their Lya 
could also be unmeasured if circular apertures are used. 

Lya emission more extended than the F UV has been 
report ed in numerous high- z samples. iFvnbo et al.l 
(|200l) remarked upon a few such objects at the brighter 
end of the luminosity distribution of the 27 narrow- 
band selected galaxies, and the extr e mely deep spec- 
troscopic observations of iRauch et al.l (|2008h uncovered 
28 Lya galaxies, ten of which were classified as ex- 
tended. Samples of Lya blobs (e.g. iMatsuda et al.l 120121 : 
iPrescott et~aLll2012[ ) may be many times the size of their 
counterpart galaxies, if indeed counterparts are identified 
at all. Here we report that every galaxy in the sample 
that emits Lya does so by producing a halo; on average 
the halo is over twice the linear size of Ha and the FUV. 

By stacking nar r owban d images of LBGs at (z) = 
2.65, iSteidel et al.l (|2011l) reported Lya halos that ex- 
tend many tens of kpc, probably probing the neutral 
circumgalactic medium (CGM) out to the virial radius. 
Subdividing the full sample by Lya properties, the halos 
at radii larger than 20-30 physical kpc show very sim- 
ilar scale lengths in all subsamples (although different 
central surface brightnesses) , even when central Lya ab- 
sorption is found. At small radii the subsamples exhibit 
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Table 1 

The LARS sample: properties and sizes. 



LARS 


Common name 


R.A. 


Dec. 


z 


ipuv 


pFUV 
■ rl P20 


-"-P20 


r)Lyct 
n P20 


€p< 2 ) 




/3-slope 


Ha/H/3 


ID 




h:m:s 


d:m:s 




Lq 


kpc 


kpc 


kpc 


kpc 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


01 


Mrk259 


13:28:44.0 


+43:55:49.9 


0.028 


9.92 


1.18 


1.29 


4.36 


7.87 


3.37 


-1.83 


3.08 


02 




09:07:04.9 


+53:26:56.5 


0.030 


9.48 


1.12 


1.17 


2.67 


8.41 


2.27 


-2.02 


3.08 


03 


Arp238 


13:15:35.1 


+62:07:27.2 


0.031 


9.52 


0.84 


0.97 


0.75 


8.68 


0.77 


-0.57 


5.18 


01 




13:07:28.2 


+54:26:50.7 


0.033 


9.93 


3.79 


1.57 




9.22 




-1.76 


3.48 


05 


Mrk 1486 


13:59:51.0 


+57:26:23.0 


0.034 


10.0 


0.93 


1.24 


3.24 


9.49 


2.61 


-2.09 


3.06 


06 


KISSR2019 


15:45:44.5 


+44:15:49.9 


0.034 


9.20 


3.65 


0.66 




9.48 




-1.85 


2.96 


07 


IRAS 1313+2938 


13:16:03.9 


+29:22:54.2 


0.038 


9.75 


0.85 


0.89 


3.01 


10.5 


3.37 


-1.94 


3.37 


08 




12:50:13.7 


+07:34:44.2 


0.038 


10.2 


5.01 


3.89 


4.35 


10.5 


1.12 


-0.90 


4.09 


09 


IRAS 0820+2816 


08:23:54.9 


+28:06:22.8 


0.047 


10.5 


5.00 


4.21 


>12.0 


12.9 


>2.85 


-1.52 


3.48 


10 


Mrk 0061 


13:01:41.5 


+29:22:53.2 


0.057 


9.74 


2.34 


2.63 


5.49 


15.5 


2.08 


-1.36 


3.93 


11 




14:03:47.1 


+06:28:15.0 


0.084 


10.7 


8.00 


6.81 


15.5 


22.1 


2.27 


-1.50 


4.60 


12 


SBS 0934+547 


09:38:13.5 


+54:28:25.3 


0.102 


10.5 


1.78 


2.03 


7.06 


26.3 


3.48 


-1.92 


3.21 


13 


IRAS 0147+1254 


01:50:28.4 


+ 13:08:59.2 


0.147 


10.6 


3.83 


4.68 


8.12 


36.0 


1.74 


-1.53 


4.07 


11 




09:26:00.3 


+44:27:36.0 


0.181 


10.7 


0.79 


1.62 


5.86 


42.7 


3.62 


-2.22 


3.13 



Note. — Coordinates (3 and 4) arc J2000. Rcdshifts (5) are derived from SDSS. (6) arc log(i/Ljy) in solar luminosities. (7—9) arc the Pctrosian 
radii with rj — 0.2, Rp20- (10) shows the physical scale corresponding to an angular size of 14 arcsec at the redshift of each galaxy - this corresponds 
to half the diametric size of the ACS/SBC (28 arcsec / 2) and describes the maximum usable scale to which we can probe Lya. (11) gives the 
relative extension of Lya relative toHa. (12) gives the UV slope, /3, derived from HST imaging. (13) gives the Ha/H/3 ratio, derived from SDSS 
spectroscopy. 
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Figure 3. Comparison of the Petrosian radii (rj = 0.2), i?p2Cb in continuum-subtracted Lya, Ho, and the UV continuum. In cases where a 

galaxy is a net Lya absorber its size is undefined, and ^ptq has been set to a small negative value - it could in principle also be very large. 
The Left panel shows how Lya sizes compare with the FUV, which can be similarly derived at high-z. The Right panel makes the same 
comparison against R^20' wn i crl directly contrasts the intrinsic and emitted Lya sizes. When net Lya emission is found it is systematically 
extended, taking mean values of R^/R^=2.9 (left) and R^/R^ =2A (right). 



profiles that differ markedly, dropping rapidly to ~ for 
the Lya-absorbing sample but steepening by varying de- 
grees in all others. Even the steepest central profiles, 
however, still run much flatter than those of the stel- 
lar continuum, and this change likely marks the onset of 
higher density gaseous disks or similar. 

From the various z « 2.7 Lya profiles of iSteidel et al.l 
(|2011h we calculate i?p2o usm S the same method as for 
our sample, and dividing by -Rparj from the continuum 
profile we obtain ^Lya (now relative to the UV). These 
raw values range between £Lya=3-8 for the non-LAEs, 
and 5.9 for the LAE-only sample, and are notably big- 
ger than our largest ^Lya- However, under the assump- 
tion that the inner and outer profiles mark physically 
different regimes that may not be the same in low-z 
galaxies, we also s ubtract the exponential halo fits of 
ISteidel et al.l ()2011h and repeat the exercise; this yields 
a range of ^Lya=0.84 to 2.0. This is now smaller than 
many of our values, although close to the average and 



the dispersion of the high-z sample is obviously lost in 
the stacking proc ess. On the o t her h and, the UV con- 
tinuum profile of ISteidel et al.l (| 2 1 If) is dominated by 
atmospheric seeing. If we instead use the continuum ef- 
fective radius of BM/BX galaxies and LBGs from HST 
imaging (IMosleh et alJl201lh we compute -Rp2o V ~ ^ kpc, 
which would increase all the ^Lya c|uoted above by a fac- 
tor of 2.5. ^Lya from the raw data would then become 
much larger than we measure in the local universe (up to 
15), and ^Lya in halo-subtracted profiles that are roughly 
consistent (2.1 to 4.8). 

LARS observations probe scales far below the tens of 
kpc sampled at high-z on a case-by-case basis. The 
galaxies likely include the range between , or roughly 
bracke ting, the averaged subsamplcs of ISteidel et al.l 
(| 2 1 If) . Our imaging also suggests this extension to be a 
very common property of Lya-cmitting galaxies, and its 
onset begins almost immediately in the inner few kpc: 
we find seven galaxies with FUV Petrosian radii below 
2 kpc, five of which have corresponding Lya radii three 
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Figure 4. Correlations between the Lya vs. Ha extension, ^L ya , and measures of the galaxy reddening. The Left panel shows the UV 
slope, measured from HST imaging, while the Right panel shows the nebular attenuation measured from SDSS spectroscopy at Ha and 
H/3. Net absorbing galaxies are set to zero and ringed, but could in principle also be very extended. Rank correlation coefficients of the 
Spearman p test and the associated probability of the no-correlation hypothesis (accounting for ties and the small sample) arc given in the 
top-right corners. Symbols are the same as Figure [3] 

Ha /H/3 ratio is the canonical probe of nebular reddening 
(i.e. that which is to zeroth order expected for Lya) used 
in studies of low-z and Galactic nebulae. /3 and Ha /H/3 
are listed in Table [TJ 

Both measures of dust content strongly anti-correlate 
with ^Lya although the sample is small (N — 12 de- 
fined sizes in Lya). To assess its significance we com- 
pute the Spearman rank correlation coefficient, p, which 
yields p = —0.73 and —0.61 for the anti-correlation of 
^Lya with /3 and Ha/H/3, respectively. This corresponds 
to likelihoods of the null hypothesis - that this correla- 
tion arises purely by chance - amounting to 0.7 percent 
(UV slope), and 3.6 percent (Ha/H/3). 

The halo-dust phenomenon appears not to be a direct 
effect of radiative transfer. We have performed new tes t 
simulations with the McLya code (IVerhamme et al.l2006f) . 
by tuning the gas-t o -dust ratio in the synthetic galaxy 
of IVerhamme et al.l (|2012l ) . Indeed the surface bright- 
ness does scale with dust abundance but the light profile 



times larger. 

It is also noteworthy that iSteidel et al.l (| 2 1 lh find dif- 
ferent median dust attenuations for the Lya-emitting 
a nd non-emitting sub samples, almost precisely as we did 
in lHaves et al.l (|2010l ). LAEs, which show extended cen- 
tral peaks, were determine d to have stellar E b-v^0-09 
magnitudes (c.f. 0.085 in lHaves et al.l 120101) while ab- 
sorbers show .Eb_v=0.19 (c.f. 0.23 for ou r Ha-selected 
sampl e) . Adopting the prescription of IMeurer et al.l 
(|1999| ) the stellar Eb-v measurements for the Steidel 
et al. samples correspond to j3 sloped of —1.77 (LAEs) 
and —1.27 (Lya absorbers). Bluntly accounting for a 
factor of 2.27 that connects st ellar Eb-v to its ne bular 
equivalent in local starbursts (|Calzetti et al.l [20001) . the 
same stellar Eb-v would equate to Ha/H/3 ratios of 3.5 
(LAEs) and 4.4 (absorbers). In the next Section we will 
show case-by-case that Lya halos systematically become 
more extended with decreasing dust contents. 

5. LYMAN ALPHA EXTENSION AND DUST CONTENTS 

In Hayes et al. (in prep) we compute many global 
properties for the sample, in order to study the processes 
complicit in Lya transport. Indeed that paper will in- 
clude a complete analysis of correlations between Lya 
transmission, halo sizes, and many other properties; in 
this Letter we restrict ourselves to observables that scale 
with the dust content. It is noteworthy for the moment, 
however, that we find no correlation between ^Lya and 
the total Lya luminosity. In Figure [4] we show how £,L ya 
compares with both the UV continuum slope /3 and the 
Ha/H/3 ratio. We note that the SDSS fibers are on aver- 
age smaller than the Lya radii, but do capture the bulk 
of the nebular emission, and fluxes can easily be mea- 
sured with out contamina t ion of [Nil] and stellar absorp- 
tion. Since IMeurer et al.l (|1999l ) /3 has been used almost 
ubiquitously as a proxy of stellar attenuation in high-z 
galaxies; here we measure /3 from aperture-matched HST 
imaging using the FUV (SBC/F140LP or F150LP) and 
the [/-band (UVIS/F336W or F390W) filters. With col- 
ors between /3 m —2.2 and —0.6 our objects have similar 
UV slopes to the vast majority of those fo und in z = 2 — 4 
Lya-emitting galaxies dBlanc et al.ll20lil ). Similarly the 

16 UV continuum flux density, parameterized by a power-law of 
the form fx oc \P . 



(therefore -Rp2o ) does not, and the ^Lya _ dust trend must 
be a secondary correlation. A scenario is needed in which 
galaxies decrease the relative size of their Hi envelopes 
as the absolute dust content increases. A sequence in 
which neutral gas settles into the galaxy (reducing ^Lya) 
and subsequently forms stars (creating more dust) would 
explain the trend, but without yet having obtained spa- 
tially resolved Hi data this is conjecture. 

Scattering also has the potential to spread Lya over 
such an area that its surface brightness decreases greatly. 
In such a case, scattered radiation measured at large radii 
may not be sufficient to recover flux from a broad central 
absorption, making ^Lya observationally undefined when 
it is actually very large. The trend of ^Lya increasing 
in bluer galaxies, then, is also able to explain the unde- 
fined sizes of LARS 04 and 06, at their measured dust 
abundance. Similar considerations would also explain 
the non-detection of Lya in local gas-rich but metal- 
and dust-poor dwarf starbursts such as I Zw 18 and 
SBS 0335-052 (IKunth et al.H1994tlMas-Hesse et al.ll200l 
IQstlin et aLll2009D . as discussed in lAtek et all (I2009bfl. 

We have empir ically shown before (jAtek et alj|2009at 
lHaves et al.ll2010T ) that the global escape fraction of Lya 
photons anti-corre lates strongly with attenuation (also 
IKornei et al1[20Tfj| in LBGs). We now demonstrate that 
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at lower Eb-Vi the more strongly emitting galaxies are 
likely to also spread their Lya over larger surfaces. Thus 
while they do transmit more of their Lya, it may be that 
more of the transferred Lya is observationally lost out- 
side photometric apertures. This may also explain the 
lac k of correlation b e tween Lya/H/3 and Eb-v observed 
bv iGiavalisco et al.l (|1996l) . compared to trends seen in 
other samples: the aperture of the IUE probed just 3 kpc 
at z = 0.01 and if more Lya is lost in bluer galaxies 
the Lya/Balmer ratios would be artificially lowered in 
in such systems. This could in part mask an underlying 
correlation. By a similar token, galaxies that can very 
efficiently scatter Lya photons may not be recovered at 
all, despite frequently showing very blue UV colors. De- 
termining precisely how Lya profiles are modified for a 
given set of host properties will provide a cornerstone for 
interpreting future large high-z surveys. 
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